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SUMMARY 

The  NAE  Twin  Otter  atmospheric  research  aircraft  carries 
instruments  mounted  under  the  wings  to  count  and  image  cloud  and  precipi¬ 
tation  particles  for  cloud  physics  studies.  As  part  of  an  investigation  on  the 
influence  of  the  aircraft  presence  on  these  measurements,  airflow  velocities 
about  the  wing  have  been  calculated.  These  theoretical  results  have  been 
compared  with  in-flight  measurements  taken  by  pressure  probes  mounted  in 
a  cannister  suspended  from  a  long  and  a  short  pylon  below  the  wing  at  two 
spanwise  locations. 

The  Twin  Otter  wing  has  been  represented  by  a  Joukowski  aero¬ 
foil  of  the  same  maximum  thickness  but  with  a  cusped  trailing  edge.  The 
wing  camber  has  been  neglected  but  aspect  ratio  and  planform  effects  have 
been  included.  The  effect  of  the  cannister  was  included  in  the  comparison  of 
theory  with  experiment.  Use  of  the  theory  aft  of  the  quarter  chord  will 
require  experimental  verification. 

Both  experimental  and  theoretical  results  for  the  local  flow 
velocity  ratio  are  linear  functions  of  aircraft  lift  coefficient  (CL )  in  the 
range  of  CL  between  0.35  to  1.06.  The  worst  discrepancy  between  theory 
and  experiment  in  the  CL  range  of  0.4  to  0.6  is  -2.8%  to  +1.1%.  The  results 
for  the  inboard /long  pylon  and  outboard /short  pylon  have  considerably  less 
deviation. 


The  theory  has  been  experimentally  verified  to  the  above  accuracy 
and  can  be  used  to  describe  the  spatial  variation  of  the  flow  velocity  ahead 
of  the  wing.  Applications  of  the  theory  to  the  study  of  aircraft  flow  effects 
on  cloud  physics  measurements  can  proceed  with  a  quantifiable  level  of 
confidence. 


I  Accession  F*r 


L’avion  Twin  Otter  pour  l’etude  de  l’atmosphere  de  TEtablissement 
aeronautique  national  transporte  sous  ses  ailes  des  instruments  qui  comptent 
et  represented  par  une  image  les  particules  de  nuage  et  de  precipitation  pour 
l’etude  de  la  physique  des  nuages.  Dans  le  cadre  d’une  recherche  sur  l’influ- 
ence  de  l’avion  sur  ces  mesures,  on  a  calcule  les  vitesses  de  Pecoulement 
aerodynamique  sur  les  ailes.  Ces  resultats  theoriques  ont  ete  compares  a  des 
mesures  prises  en  vol  a  l’aide  de  sondes  de  pression  montees  sur  une  boite 
suspendue  par  un  mat  long  et  par  un  mat  court  en  deux  endroits  dans  le  sens 
de  l’envergure  de  l’aile. 

L’aile  du  Twin  Otter  a  ete  representee  par  un  profil  aerodynamique 
de  Joukowski  de  la  meme  epaisseur  maximale,  mais  comportant  un  bord  de 
fuite  redente.  On  n’a  pas  tenu  compte  de  la  courbure  de  l’aile,  mais  on  a 
inclus  les  donnees  des  effets  de  l’allongement  et  de  la  forme  en  plan.  On  a 
tenu  compte  de  l’effet  produit  par  la  boite  dans  la  comparaison  entre  la 
theorie  et  l’experimentation.  L’application  des  donnees  theoriques  a  Tarriere 
du  quart  avant  doit  etre  effectuee  apres  une  comparaison  plus  approfondie 
avec  les  donnees  de  Pexperimentation. 

Les  resultats  theoriques  et  ceux  qui  decoulent  de  Texperimenta- 
tion,  en  ce  qui  a  trait  au  rapport  de  la  vitesse  d’ecoulement  local,  sont  des 
fonctions  lineaires  du  coefficient  de  portance  (Cz)  dans  la  plage  de  cette 
demiere  comprise  entre  0.35  et  1.06.  Le  pire  ecart  entre  la  theorie  et  l’expe- 
rimentation  dans  la  plage  de  Cz  comprise  entre  0.4  et  0.6  varie  de  -2.8%  a 
+  1.1%.  Les  resultats  obtenus  pour  le  mat  interieur  long  et  le  mat  exterieur 
court  presentent  un  ecart  beaucoup  moins  important. 

La  theorie  a  ete  verifiee  par  I’experimentation  au  niveau  de  preci¬ 
sion  indique  ci-dessus,  et  elle  peut  etre  employee  pour  decrire  la  fluctuation 
spatiale  de  la  vitesse  d’ecoulement  a  l’avant  de  l’aile.  On  peut  appliquer  la 
theorie  a  l’etude  des  effets  de  l’ecoulement  sur  les  mesures  en  physique  des 
nuages  avec  une  certaine  confiance. 
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THEORETICAL  AND  MEASURED  AIRFI  OW  AHOlJT  THE  TWIN  OTTER  WING 


INTRODUCTION 

Joukowski  aerofoil  theory  is  applied  to  study  the  spatial  variation  of  the  velocity  field  in 
the  vicinity  of  the  leading-edge  of  the  Twin  Otter  wing.  The  NAE  Twin  Otter  atmospheric  research 
aircraft  is  instrumented  to  measure  and  to  compute  in  real  time  atmospheric  motion  and  state 
parameters.  For  several  years  it  has  been  used  in  a  cloud  physics  research  role  in  co-operation  with  the 
Atmospheric  Environment  Service.  Particle  measuring  probes  are  carried  under  the  aircraft  wings  to 
measure  cloud  droplet  spectra,  ice  particle  shapes,  etc.  These  probes  are  mounted  a  short  distance 
ahead  of  cannisters  which  are  suspended  by  pylons  from  the  wing  lower  surface.  There  is  evidence  to 
suggest  that  images  and  measured  concentrations  of  cloud  droplets  and  ice  particles  may  be  affected 
by  the  flow  around  the  wings.  It  was  this  concern  that  motivated  the  initiation  of  the  work  described 
here. 

The  scope  of  the  report  is  limited  to  the  establishment  of  a  model  for  predicting  the  local 
flow  velocity  near  the  wing  as  a  function  of  position.  Applications  of  the  model  to  cloud  physics  will 
be  the  subject  of  another  report. 

No  new  theory  will  be  developed  here.  The  results  of  Glauert  (Ref.  1)  and  Abbot  and 
von  Doenhoff  (Ref.  2)  are  used  as  a  theoretical  basis.  A  summary  of  the  theory  is  given  in  Appendix  I. 
Kuethe  and  Schetzer  (Ref.  3)  have  a  very  clear  development  of  some  of  the  necessary  concepts  of 
wing  theory  that  are  required  in  order  to  apply  the  two-dimensional  results  of  Glauert. 

Glauert  (Ref.  1)  first  published  his  theory  in  book  form  in  1926  and  a  slightly  extended 
version  was  published  in  1947.  He  draws  upon  theory  developed  in  the  period  1902-1924.  We  will 
refer  to  his  results  and  quote  appropriate  page  numbers  if  necessary.  Abbot  and  von  Doenhoff  (Ref.  2) 
also  rely  heavily  on  the  theory  of  Glauert  but  their  treatment  of  the  subject  has  exceptional  clarity  of 
exposition.  Their  work  was  first  published  in  1949;  the  1958  edition  quoted  here  is  identical  except 
for  a  few  corrections. 

It  is  clear  that  the  theory  being  applied  here  is  not  new.  In  this  work  the  full  generality  of 
the  method  is  not  utilized.  The  simplest  transformation  formula  is  used  and  wing  camber  is  neglected. 
The  trailing  edge  of  the  wing  is  allowed  to  be  a  cusp  and  no  attempt  has  been  made  to  include  the 
correct  trailing  edge  angle.  This  latter  restriction  is  thought  to  be  not  serious  because  only  flow  ahead 
of  the  wing  quarter -chord  point  is  of  interest  in  the  present  application. 

One  of  the  main  contributions  of  this  work  is  the  experimental  verification  of  the  method 
at  the  locations  near  the  wing  where  the  particle  measuring  probes  are  mounted.  Without  this  con¬ 
firmation,  application  of  the  theory  would  be  of  academic  interest  only.  Another  important  result  is 
the  linear  correlation  of  the  magnitude  of  the  local  flow  velocity  with  aircraft  lift  coefficient  in  both 
theoretical  and  experimental  instances. 

The  next  section  discusses  application  of  the  conformal  transformation  described  in 
Appendix  I  to  the  Twin  Otter  wing.  Some  minor  extensions  to  Glauert ’s  results  are  included. 

GEOMETRICAL  CONSIDERATIONS 

Using  Equation  (1-5)  from  the  summary  of  theory  in  Appendix  I,  the  trailing  edge  of  the 
wing  transforms  from  ‘a’  in  the  z  plane  to  2a  in  the  f  plane.  In  the  z  plane,  the  leading  edge  is  located 
at  -2(a+e)  relative  to  the  trailing  edge.  Hence,  the  z  co-ordinate  of  the  leading  edge  is  -(a+2e).  Using 
Equation  (1-5),  the  leading  edge  in  the  wing  plane  is  located  at 


-3- 


The  Twin  Otter  wing  has  an  NACA  0016  thickness  distribution  which  is  represented 
analytically  by  the  following  equation  (Ref.  2,  p.  113): 

rjl  =  5(t/c)  |o.2969£*  -  0.126  ^  -  .3516  +  .2843  -  .1015  (4) 

where  rjj  and  are  fractions  of  the  chord  c.  The  thickness  at  any  chordwise  location  is  obtained  by 
doubling  the  result  of  (4).  Note  that  with  =  1.  the  trailing  edge  thickness  is  not  zero.  Figure  1  shows 
the  Twin  Otter  wing  shape  obtained  by  Equation  (4)  and  the  Joukowski  aerofoil  profile  from  Equa¬ 
tions  (2)  and  (3).  The  probe  locations  and  the  Joukowski  parameter  values  are  also  shown.  There  is 
extremely  good  agreement  between  the  two  shapes  as  far  back  as  the  quarter  chord  point.  The 
Joukowski  profile  is  slightly  fatter  than  the  NACA  section  until  the  quarter  chord  point  whereupon 
the  NACA  section  becomes  noticeably  thicker.  The  trailing  edge  cusp  of  the  Joukowski  profile  is 
quite  visible. 

The  stagnation  point  at  the  leading  edge  is  located  by  forcing  the  terms  inside  the  square 
brackets  in  Equation  (1-8)  to  be  zero  when  the  value  of  T  from  Equation  (1-9)  is  inserted.  In  the  z 
plane,  after  some  lengthy  algebra,  the  result  is: 


=  -  ^e+(a+e)cos2oQ^  -  i(a+e)sin2c*o 


The  stagnation  point  in  the  wing  plane  is  obtained  by  substituting  the  result  of  Equation  (5)  in 
Equation  (1-5).  This  procedure  is  best  performed  numerically.  If  the  approximate  forms  for  the  aero¬ 
foil  co-ordinates  due  to  Glauert  are  used,  one  finds  that  the  stagnation  point  does  not  lie  on  the  wing. 
In  the  present  case,  it  does. 

It  is  apparent  that  the  flow  velocity  at  a  given  point  depends  only  on  aQ  and  e.  The  calcula¬ 
tion  of  the  thickness  parameter  e  was  discussed  previously.  We  now  attempt  to  obtain  a  value  of 
under  the  assumptions  that  the  flight  speed  is  not  too  fast  (M  <  .7)  and  that  the  spanwise  lift  distribu¬ 
tion  is  elliptical.  For  the  Twin  Otter,  the  Mach  number  restriction  is  easily  met  and  the  lift  distribution 
assumption  is  approximately  fulfilled  (Ref.  1,  p.  155).  The  lift  at  the  outboard  probe  location  will  be 
slightly  under-estimated  but  there  is  no  significant  discrepancy  at  the  inboard  location. 

The  real  wing  is  describable  only  by  three-dimensional  flow  while  the  potential  function 
(Eq.  (1-8))  is  defined  for  two-dimensional  flow.  Kuethe  and  Schetzer  (Ref.  3,  p.  96)  utilize  two- 
dimensional  relations  to  obtain  three-dimensional  results  for  the  wing:  “The  analysis  proceeds  on  the 
assumption  that  the  flow  about  each  unit  span  of  the  finite  wing  is  two-dimensional  .  .  .”.  We  will 
use  the  converse  statement  and  obtain  two-dimensional  values  from  three-dimensional  wing  results. 

The  equation  for  the  three-dimensional  lift  curve  slope  (m)  in  terms  of  the  two-dimensional  result 
(mQ )  (Ref.  3,  p.  100)  is: 


1  + -  (1+r) 

ttAR 


where  AR  is  the  wing  aspect  ratio  and  r  accounts  for  planform  effects.  Values  of  r  are  tabulated 
(Ref.  1,  p.  147)  for  various  values  of  AR/m  and  taper  ratios. 


The  section  lift  coefficient  Ce  is  defined  as 


c.  =  e/(*pv2c) 


This  can  be  related  to  the  sectional  circulation  by  the  Kutta-Joukowski  Theorem  (Ref.  2,  p.  67) 
which  results  after  integration : 


s  =  pvr  (f 

where  the  value  of  T  is  obtained  from  Equation  (1-9).  Performing  the  necessary  substitutions  and 
simplifying  leads  to : 


Cg  =  2  7T  (1+e/a)  sin  aQ  (! 

The  chord  has  been  taken  to  be  4a.  Differentiating  Equation  (9)  with  respect  to  oo  results  in  an 
expression  for  the  two-dimensional  lift  curve  slope  for  values  of  cxq  small  enough  to  justify  the  ap¬ 
proximation  cos  aQ  3  1.  The  result  is: 


mQ  =  2  7r  (1+e/a)  (1( 

Using  the  value  of  e/a  obtained  from  Table  I  at  16%  thickness  yields  the  value  of  mo  as 
7.155  per  radian.  The  aspect  ratio  of  the  Twin  Otter  wing  is  10.  These  two  numerical  values  were 
applied  to  the  data  for  r  in  Reference  1  which  allowed  the  value  of  r  to  be  estimated  as  0.212. 

The  three-dimensional  lift  curve  slope  can  be  calculated  from  Equation  (6)  using  the 
necessary  numerical  values.  The  value  of  m  is  5.607  per  radian  or  .0979  per  degree.  Elliptic  loading 
requires  that  the  product  (a^C^c)  must  vary  elliptically  (Ref.  3,  p.  99).  Therefore,  for  a  rectangular 
wing  with  no  spanwise  thickness  variation  or  twist,  the  angle  of  attack  «o  must  be  an  elliptical  func¬ 
tion  of  the  span: 


a0(s)  =  «c  v  i-s2  (i: 

where  s  is  the  non-dimensional  spanwise  co-ordinate  and  ac  is  the  value  of  aQ  at  the  wing  center  line. 

The  aircraft  CL  is  known  for  a  given  flight  condition.  The  Prandtl-Glauert  rule  is  applied 
to  m  to  obtain  the  lift  curve  slope  m'  when  the  Mach  number  is  not  zero: 


r 

m 


m 


(12 


This  result  is  in  “remarkably  good  agreement”  (Ref.  2,  p.  256)  with  experimental  measurements  up 
to  .7. 


The  average  wing  angle  of  attack  a  is  obtained  from 


(1! 


The  value  of  c*c  is  chosen  to  give  the  same  average  angle  of  attack  as  a: 

4  _ 

a  =  -  a  (1- 

c  it 

The  local  angle  of  attack  aQ  is  now  obtained  from  Equation  (11)  at  the  appropriate  span- 
wise  location.  The  calculation  of  the  velocity  field  about  the  wing  is  now  possible. 


The  Twin  Otter  wing  has  a  small  camber  but  the  details  of  the  shape  of  the  chord  line  were 
not  readily  available.  The  effect  of  adding  a  small  amount  of  camber  to  a  symmetrical  wing  is  to 
increase  the  lift  coefficient  without  changing  the  lift  curve  slope  (Ref.  3,  see  18.5). 

In  fact  Glauert  (Ref.  1,  p.  86)  shows  that  application  of  the  Kutta  Condition  to  a  thin 
wing  with  small  camber  results  in 


rc  =  4jra  V  sin  (cco  +  0)  (15) 

where  0  =  tan'1  (2  *camber).  If  0  is  negative,  the  wing  is  said  to  be  “reflexed”.  For  the  Twin  Otter, 
the  camber  parameter  0  should  be  positive  and  small.  The  complex  potential  from  Equation  (1-6) 
should  be  modified  to  include  camber  properly . 

However,  cambered  Joukowski  aerofoils  have  circular  arcs  for  their  mean  chord  lines.  This 
shape  does  not  fit  the  Twin  Otter  mean  chord  line  very  well.  Hence,  inclusion  of  camber  into  the 
complex  potential  does  not  seem  to  be  warranted.  Therefore,  the  effect  of  camber  on  the  local  flow 
velocity  away  from  the  wing  is  neglected. 

In  summary,  the  velocity  in  the  wing  plane  is  obtained  from  the  following  equation: 


u 


W 


V 


exp(-io:o) 


(1+e/a)2 

(z/a+e/a)2 


exp(iaQ)  +  i 


2(1 +e /a) 
(z/a+e/a) 


sin(aQ ) 


(z/a)2 
_(z/a)2  -1  _ 


(16) 


Equation  (1-7)  should  be  used  to  obtain  z  for  a  desired  f .  The  magnitude  of  the  local  flow  velocity 
relative  to  the  value  at  infinity  is: 


v  ■  7  vV+»„J  an 

From  the  previous  discussion,  after  the  wing  thickness,  aspect  ratio  and  planform  have  been 
given,  the  flow  velocity  depends  only  on  z/a  and  aQ  The  latter  is  determined  from  the  aircraft  CL  and 
to  a  much  lesser  extent  from  the  Mach  number  as  well. 

For  a  symmetrical  wing  at  angles  of  attack  small  enough  for  first  term  approximations  to 
the  trigonometric  functions  to  be  acceptable,  both  uw  /V  and  vw  /V  are  linear  functions  of  ao : 

uw 

—  =  f,  (x,y)«0 


■y  -  f2  (x«y>ao 

Substitution  of  these  results  in  Equations  (16)  and  (17)  and  differentiating  with  respect 
to  a0  shows  that  9V/3cc0  is  a  constant.  The  aircraft  lift  coefficient  is  a  linear  function  of  x0  if  com¬ 
pressibility  effects  are  temporarily  ignored.  Hence,  we  would  expect  a  plot  of  V_versus  CL  to  be 
represented  very  closely  by  a  straight  line  at  low  Mach  numbers.  Computations  of  V  versus  CL  were 
performed  from  Equations  (16)  and  (17).  Straight  lines  fitted  the  data  well  for  values  of  a0  below 
the  stall. 


-  6  - 


The  spatial  variation  of  V  has  been  calculated  for  a  few  values  of  CL  with  compressibility 
effects  being  temporarily  ignored .  All  parameters  have  been  chosen  to  suit  the  Twin  Otter.  Lines  of 
constant  V  (isovels)  are  shown  in  Figures  2  to  5  for  the  inboard  probe  location.  The  location  of  the 
measuring  probes  at  each  pylon  length  under  the  wing  are  shown. 

As  CL  increases,  the  stagnation  point  moves  farther  around  the  underside  of  the  wing.  The 
isovels  are  very  nearly  symmetrical  about  the  leading  edge  at  low  CL . 

Consider  Figure  4  which  shows  the  isovels  at  CL  =  0.55.  Moving  from  left  to  right  at 
constant  17/a,  the  velocity  variation  is  markedly  different  depending  on  the  value  of  17 /a.  When  the 
latter  variable  is  -.4  (say)  initially  one  sees  V  slightly  less  than  1.  Moving  to  the  aft  limit  one 
encounters  a  decreasing  velocity  and  then  an  increasing  velocity.  The  gradients  are  extremely  high 
near  the  stagnation  point.  However,  if  17 /a  =  -1.4,  the  velocity  does  not  change  nearly  as  much  and 
the  gradients  are  much  smaller.  As  CL  increases,  velocity  decrements  are  felt  farther  away  from  the 
wing. 

The  stream  function  \p  is  the  imaginary  part  of  the  complex  potential.  It  is  quite  straight¬ 
forward  to  calculate  the  equation  for  a  streamline  (i j/  =  constant)  that  passes  through  a  given  point 
near  the  wing.  Figures  6(a)  and  6(b)  show  two  such  lines  for  high  speed  and  low  speed  flight  respec¬ 
tively  that  pass  through  the  same  point  under  the  wing.  The  scales  of  17 /a  are  different.  The  stream¬ 
lines  are  each  asymptotic  to  the  appropriate  a0 ,  but  near  the  wing,  there  is  an  upwash  ahead  of  the 
leading  edge.  The  larger  the  CL ,  the  larger  the  upwash.  The  flow  then  slopes  either  down  or  further 
up  depending  on  the  value  of  CL  .  The  streamline  to  the  stagnation  point  ( 1//  =  0)  is  shown  in  each  case. 

The  magnitude  of  the  flow  velocity  along  the  streamlines  of  Figures  6(a)  and  6(b)  are 
shown  in  Figures  7(a)  and  7(b).  For  a  low  value  of  CL ,  (Fig.  7(a))  the  flow  velocity  is  not  very  dif¬ 
ferent  from  the  flow  velocity  at  infinity.  There  is  a  local  minimum  just  ahead  of  the  leading  edge.  At 
higher  CL ,  (Fig.  7(b)),  the  minimum  is  more  pronounced  and  occurs  a  little  farther  to  the  rear.  The 
minimum  is  much  lower  and  the  local  flow  reaches  the  aircraft  speed  much  farther  upstream. 

Figures  8(a)  and  8(b)  show  the  local  flow  angle  along  the  same  streamlines  as  in  Figures  6(a) 
and  6(b).  In  the  vicinity  of  the  leading  edge,  the  flow  angles  exhibit  large  gradients.  There  is  a  local 
maximum  near  the  leading  edge  in  each  case.  The  lower  the  value  of  CL ,  the  smaller  is  the  change  in 
local  flow  angle  compared  to  aQ . 

MEASUREMENT  OF  FLOW  VELOCITY 

A  pitot/static  and  a  total  pressure  probe  were  installed  side-by-side  on  the  nose  of  a  standard 
Particle  Measuring  Systems  long  cannister.  Figures  9(a)  and  9(b)  show  this  cannister  mounted  on  a 
short  and  on  a  long  pylon  beneath  the  Twin  Otter  wing.  The  relevant  dimensions  are  given  in 
Appendix  II. 

It  was  desired  to  measure  the  static  pressure  (Ps )  and  dynamic  pressure  (PT  -Ps )  ahead  of 
the  cannister  in  the  precise  sampling  plane  of  a  cloud  particle  measuring  probe,  but  this  was  not 
perfectly  achievable  because  of  physical  limitations  in  mounting  the  available  pressure  probes.  Only 
the  static  pressure  from  the  pitot/static  probe  was  used.  The  centre  of  its  two  rows  of  pressure  ports, 
with  axial  displacement  of  V4  inch,  was  xh  inch  ahead  of  the  particle  sampling  plane,  and  the  inlet  of 
the  short  total  pressure  probe  was  %  inch  aft  of  this  plane.  The  result  is  that  the  static  pressure  meas¬ 
urements  were  made  at  a  point  l'A  inches  ahead  of  the  total  pressure  measurements  and,  on  average, 

V4  inch  ahead  of  the  sampling  plane.  Since  static  pressure  is  the  critical  parameter  (PT  should  be  con¬ 
stant),  it  is  considered  to  be  measured  within  an  acceptable  dimensional  tolerance  (0.3%  of  chord); 
hence,  its  measurements  can  be  considered  to  apply  at  the  sampling  plane.  The  2  inch  lateral  separa¬ 
tion  between  the  static  and  the  total  pressure  probes  can  be  ignored,  because  there  is  negligible  span- 
wise  variation  in  the  flow  about  the  constant  section  wing  of  the  Twin  Otter  at  the  probe  mounting 
locations. 


Static  and  dynamic  pressure  transducers  were  housed  within  the  cannister.  The  transducers 
had  very  short  pressure  lines  to  the  probes.  The  dynamic  pressure  was  measured  differentially  between 
the  total  head  probe  (PT )  and  the  static  pressure  (Ps )  from  the  pitot/static  probe.  Signals  from  the 
transducers  were  carefully  calibrated  and  recorded  digitally  at  16  samples/second  on  the  NAE  record¬ 
ing  system  along  with  all  the  other  parameters  normally  recorded  on  the  Twin  Otter  (Ref.  7).  The 
recording  resolutions  for  Ps  and  PT-PS  were  0.97  and  .07  millibars  per  bit  respectively,  while  RMS 
errors  in  the  quadratic  fits  to  the  calibrations  were  0.73  and  .02  millibars. 

Four  test  flight  configurations  were  flown.  The  airspeed  measuring  cannister  was  mounted 
on  both  short  and  long  pylons,  and  at  both  the  inboard  and  outboard  positions  on  the  wing.  Stabi¬ 
lized  level  flight  was  selected  at  several  altitudes  over  a  range  of  indicated  airspeeds  from  near  the 
flaps-up  stall  speed  to  near  full  speed  at  about  10  knot  increments.  For  each  test  run  the  event 
marker  on  the  recording  system  was  depressed  for  about  10  seconds  to  denote  a  valid  stabilized 
flight  condition.  During  the  subsequent  analysis,  data  were  averaged  over  this  period.  For  at  least 
one  altitude  on  each  flight,  runs  were  flown  in  both  an  increasing  and  decreasing  order  of  speeds  to 
ensure  that  there  were  no  measurement  errors  due  to  hysteresis.  Differing  fuel  loads  were  carried  to 
cover  the  usual  operating  gross  weight  range  of  the  aircraft.  Table  II  presents  the  test  points  flown 
over  the  following  wide  range  of  flight  parameters: 

(i)  indicated  airspeeds  from  80  to  140  knots  (41  to  72  m  s'1 ) 

(ii)  altitude  from  4000  to  9500  feet 

(iii)  air  temperatures  from  -1  to  20cC 

(iv)  true  airspeeds  from  93  to  167  knots  (48  to  86  m  s*1  ) 

(v)  aircraft  weight  from  9800  to  11400  lbs.  (maximum  gross  weight  is  11579  lbs.). 


TABLE  II 

TWIN  OTTER  FLIGHT  TESTS  WITH  AIRSPEED  MEASURING  PMS  CANNISTER 


Spanwise 

Pylon 

Altitude 

Position 

Length 

Feet 

Inboard 

Long 

4000 

6000 

Short 

6000 

7000 

9500 

Outboard 

Long 

6000 

7000 

9500 

Short 

6000 

9500 

Airspeed  Range 
kt  IAS 


Aircraft  Weight 
lb. 


The  local  flow  velocity  was  determined  by  using  the  perfect  gas  law  and  isentropic  flow 
equations  (Ref.  5).  In  particular,  the  local  Mach  number  (M)  was  obtained  from  measured  total  and 
static  pressures  by  utilizing  the  relation: 


where  y  is  the  ratio  of  specific  heats  for  air 

The  speed  of  sound  (Vs)  was  calculated  from  the  measured  static  temperature  (Ts)  by  the 

equation 

Vs  =  v^RT  (19) 

where  R  is  the  universal  gas  constant.  Humidity  effects  have  been  ignored  in  this  calculation  of  the 
speed  of  sound.  The  local  flow  velocity  Vn)  was  obtained  by  multiplying  the  Mach  number  by  the 
speed  of  sound 

Vm  =  V  M  (20) 

The  aircraft  lift  coefficient  is  defined  as 

CL  =  W/('/2PVa2/c  S)  (21) 

where  S  is  the  wing  area  (420  ft.2 ).  The  measured  value  of  Ts  and  the  perfect  gas  law  were  used  to 
obtain  the  air  density  p.  Aircraft  true  airspeed  derived  from  the  nose  boom  pitot/static  system  was 
computed  by  the  real-time  microprocessor  and  was  recorded  on  the  NAE  recorder.  Aircraft  weight 
was  carefully  deduced  for  each  test  point  by  taking  account  of  fuel  used  and  aircraft  weight  at  take¬ 
off.  The  definition  in  Equation  (21)  requires  steady  level  flight;  however,  some  data  were  taken  in 
light  turbulence  which  may  account  for  some  of  the  scatter  in  the  experimental  values. 

The  total  pressure  from  the  nose  boom  was  computed  from  the  sum  of  the  recorded  static 
and  dynamic  pressures.  The  total  pressure  from  the  wing-mounted  probe  was  consistently  1  to  2  mb 
higher  than  the  nose  measurement.  This  small  systematic  error  is  probably  due  to  an  error  in  the 
calibration  of  the  nose  boom  static  pressure  transducer;  however,  the  resulting  errors  in  Mach  number 
and  aircraft  true  airspeed  are  negligible. 

The  variable  V(  is  defined : 


(  7-1  , 

l  2 


b-i) 


(1.4). 


(18) 


v,  ■  vJv,i.  <22' 

It  was  calculated  for  each  test  point  as  the  ratio  of  the  velocity  at  the  probe  to  the  aircraft  true  air¬ 
speed.  The  lift  coefficient  was  calculated  from  Equation  (21). 

The  results  are  plotted  in  the  form  of  V,  versus  aircraft  C(  in  Figures  10  to  13  for  the 
two  pylon  lengths  and  the  two  spanwise  cannister  locations.  A  least-squares  procedure  was  used  to 
fit  a  straight  line  to  the  measured  data  and  the  line  is  plotted  (dashed  line)  on  the  Figures.  The 
measured  velocity  ratios  were  well  described  by  the  straight  line  in  each  case.  Table  111  shows  maxi¬ 
mum  deviations  of  measured  points  from  the  straight  line. 


TABLE  III 


DEVIATIONS  FROM  BEST  LINE  (EXPERIMENT) 


Figure 

Span 

Pylon 

Maximum 

Deviation  7r 

10 

Inboard 

Long 

±  3/4 

11 

Inboard 

Short 

±1/2  (except  near  CL  =  1) 

12 

Outboard 

Long 

±  1/2 

13 

Outboard 

Short 

+  1-1/8, -7/8 

Symbols  on  the  graphs  correspond  to  data  taken  at  different  nominal  altitudes  and  aircraft 
weights.  There  is  clearly  no  orderly  variation  of  the  results  with  altitude.  Hence,  no  further  identifica¬ 
tion  of  the  symbols  was  thought  to  be  necessary. 


Cannister  Effects 

Normally,  the  flow  about  an  aggregate  of  bodies  is  obtained  by  summing  the  complex 
potentials  of  all  the  bodies.  If  one  body  (the  cannister)  is  very  much  smaller  than  the  other  (the  wing), 
then  velocity  gradients  due  to  flow  about  the  larger  body  may  be  extremely  small  on  length  scales 
of  the  smaller  body.  The  cannister  can  be  considered  to  be  in  a  uniform  flow  whose  magnitude  and 
direction  depend  only  on  the  position  of  the  cannister  with  respect  to  the  wing. 

The  velocity  measuring  pressure  probes  are  mounted  on  long  cylindrical  cannisters  with 
hemispherical  ends.  The  length  to  diameter  ratio  is  about  9  which  allows  the  application  of  slender 
body  theory  to  the  discussion  of  flow  effects  of  the  cannister.  This  is  necessary  because  the  measure¬ 
ments  were  taken  in  a  region  where  the  flow  is  affected  by  the  cannister. 


The  local  flow  at  the  nose  of  the  cannister  is  decomposed  into  an  axial  component  (Vcosa0  ) 
and  a  cross-flow  component  ( Vsin  ).  The  latter  causes  viscous  separation  on  the  probe  body  but  the 

axial  component  is  assumed  to  be  unaffected.  Potential  flow  is  assumed  for  the  axial  flow.  It  is  suffi¬ 
cient  to  allow  the  approximation  cos  aQ  %  1  in  the  present  case  because  the  local  flow  angles  are  small 


The  dominant  effect  of  the  cannister  is  due  to  the  blunt  nose,  which  we  shall  model  as  a 
sphere.  There  will  be  some  alleviation  of  velocity  decrements  near  the  probe  nose  because  of  the 
cylindrical  after-body  but  these  will  be  neglected.  The  velocity  field  near  a  sphere  of  radius  b  (Ref.  6, 
p.  443)  looking  directly  along  the  velocity  vector  V  is: 


u  = 


(23) 


where  n  is  the  distance  from  the  sphere  center  to  the  point  of  interest. 


Velocity  decrements  at  the  measurement  point  are  easily  calculated  by  noting  that  the 
radius  of  the  sphere  is  3.25  inches  and  that  the  measurements  of  static  pressure  are  taken  at  a  point 
6-5/16  inches  in  front  of  the  nose  of  the  cannister.  At  that  point. 


u  =  .9607V  (24) 

which  means  that_V  from  the  wing  calculations  must  be  multiplied  by  .9607  to  be  comparable  to  the' 
measured  values  Vf  .  This  procedure  was  employed  for  the  following  comparison  of  theoretical  and 
experimental  velocity  ratios. 
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COMPARISON  OF  THEORY  AND  EXPERIMENT 

At  the  measured  ,  the  theoretical  velocity  ratio  was  calculated  from  Equation  1 1(5)  using 
the  previously  noted  procedure  to  obtain  ao  and  it  was  corrected  for  cannister  effects  by  Equu- 
tion  (24).  A  straight  line  of  V  versus  was  obtained  in  the  C(  range  from  0.35  to  1  .Ob.  The  theo¬ 
retical  result  corrected  for  cannister  effects  is  shown  on  Figures  10  to  13  as  a  dotted  line.  Values  of 
V  at  CL  of  0.3,  0.6  and  0.9  were  calculated.  In  all  cases,  a  straight  line  could  be  drawn  through  the 
calculated  points.  The  least  squares  fit  to  the  experimental  data  and  each  measured  point  are  shown 
on  the  same  figures  for  comparison. 

Table  IV  shows  the  difference  between  the  intercepts  for  the  theoretical  anti  experimental 
straight  lines. 

TABLE  IV 


INTERCEPT  DIFFERENCE 


Span 

Pylon 

Theory -Experiment 

Outboard 

Long 

.0073 

Short 

.0134 

Inboard 

Long 

.0031 

Short 

.0071 

1 

Table  V  shows  the  same  comparison  for  the  slopes  of  the  theoretical  and  experimental 
straight  lines: 

TABLE  V 


SLOPE  COMPARISON 


Span 

Pylon 

Theory 

Experiment 

Outboard 

Long 

-.2963 

-.3004 

Short 

-.3640 

-.3461 

Inboard 

Long 

-.3464 

-.3264 

Short 

-.4305 

-.3863 

The  intercept  of  the  theoretical  line  is  very  slightly  larger  than  that  of  the  experimental 
line  and  the  magnitude  of  the  slope  of  the  theoretical  line  is  also  slightly  larger  in  3  out  of  the  4  cases. 
Cloud  penetration  by  the  aircraft  typically  occurs  at  an  aircraft  C,  between  0.4  and  0.6  for  the  Twin 
Otter.  In  all  cases,  the  agreement  between  theory  and  experiment  is  good  in  this  range.  The  results 
for  the  inboard  pylon  location  are  under-estimated  by  the  theory  while  the  measurements  for  the 
outboard  location  are  over-estimated  by  the  theory. 

It  is  not  possible  to  isolate  reasons  for  disagreement  between  theory  and  experiment.  Even 
though  the  experimental  results  are  subject  to  error,  it  is  fair  to  attribute  discrepancies  to  the  theo¬ 
retical  model.  Some  physical  variables  were  ignored  (trailing  edge  angle,  camber  and  aircraft  thrust 
effects  on  C(  ).  However  in  the  Cj  range  of  0.4  to  0.6,  the  experimental  results  are  well-predicted  by 
the  theory  with  the  accuracy  shown  in  Table  VI,  where 


DEVIATION  OF  THEORY  FROM  EXPERIMENT 


Span 

Pylon 

Ct. 

CO  9? 

Inboard 

Long 

.4 

-  .52 

.6 

-1.02 

Inboard 

Short 

.4 

-1.67 

.6 

-2.82 

Outboard 

Long 

.4 

+  .97 

.6 

+1.10 

Outboard 

Short 

.4 

+  .65 

.6 

+  .30 

These  data  in  Table  VI  show  that  the  theory  as  presented  here  has  been  verified  with 
sufficient  accuracy.  Its  application  can  be  extended  to  the  study  of  spatial  variation  of  local  flow 
velocity  in  regions  ahead  of  the  points  at  which  measurements  were  made.  The  theory  can  be  used 
in  choosing  the  best  location  to  mount  instrumentation  that  is  affected  by  airflow. 


CONCLUSIONS 

The  Twin  Otter  wing  has  been  represented  by  a  Joukowski  aerofoil  of  the  same  maximum 
thickness  but  with  a  cusped  trailing  edge.  The  wing  camber  has  been  neglected  but  aspect  ratio  and 
planform  effects  have  been  included.  The  effect  of  the  cannister  was  included  in  the  comparison  of 
theory  and  experiment.  Application  of  the  theory  aft  of  the  quarter  chord  will  require  more  experi¬ 
mental  data. 

Both  experimental  and  theoretical  results  for  the  local  flow  velocity  ratio  are  linear  func¬ 
tions  of  aircraft  Ct  in  the  range  of  CL  between  0.35  to  1.06.  The  worst  discrepancy  between  theory 
and  experiment  in  the  CL  range  of  0.4  to  0.6  is  -2.89?  to  +1.19?.  The  results  for  the  inboard /long 
pylon  and  outboard /short  pylon  have  considerably  less  deviation. 

The  theory  has  been  experimentally  verified  to  the  above  accuracy  and  can  be  used  to 
describe  the  spatial  variation  of  the  flow  velocity  ahead  of  the  wing.  Applications  of  the  theory  to 
the  study  of  aircraft  flow  effects  on  cloud  physics  measurements  can  proceed  with  a  quantifiable 
level  of  confidence. 
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APPENDIX  I 


A  SUMMARY  OF  THEORY 


It  is  not  the  purpose  of  this  Appendix  to  re-develop  the  results  of  Glauert.  However,  the 
inclusion  of  a  few  descriptive  paragraphs  to  introduce  the  concepts  is  thought  to  have  merit. 

A  conformal  transformation  maps  a  region  of  one  plane  into  another  plane  in  such  a  manner 
that  a  small  closed  figure  in  one  plane  is  mapped  into  a  similar  closed  figure  in  another  plane  with  a 
rotation  and  magnification  (Ref.  4,  p.  54).  The  co-ordinates  in  the  wing  plane  (f  plane)  and  in  the 
transformed  plane  (z  plane)  are  related  by  the  transformation 

?  =  h(z)  (1-1) 

The  complex  potential  w  can  be  defined  in  the  z  plane  as: 


w  =  <p  +  i  \p  . 


(1-2) 


This  function  is  shown  (Ref.  2,  Section  3.3)  to  have  its  real  and  imaginary  parts  satisfy  Laplace’s 
Equation  so  that  w  may  be  interpreted  as  describing  a  physically  possible,  irrotational,  two- 
dimensional,  fluid  flow. 

The  velocity  in  the  z  plane  is 


dw 

u  -  iv  == - 

dz 


d-3) 


and  in  the  f  plane,  the  velocity  is  dw/df.  In  terms  of  the  flow  in  the  z  plane,  the  velocity  in  the  f  plane 
is 


dz 

uw  "  ivw  =  (u  -  iv)  *  ~^T  (I'4) 

where  dz/df  is  obtained  from  the  transformation  given  in  Equation  (1-1).  This  point  is  well  described 
by  Milne-Thompson  (Ref.  6,  p.  66). 

Normally,  the  flow  in  the  z  plane  is  known  and  if  one  can  find  a  suitable  transformation 
h(z),  then  the  flow  in  the  real  (f)  plane  can  be  determined. 

The  flow  around  a  circular  cylinder  of  radius  a  with  circulation  1'  is  (Ref.  2,  p.  49) 

a2  \  ir 

+  —  )  +  —  t!n(z/a) 
z  /  2rr 

The  circulation  is  related  to  the  lift  produced  by  a  flow  of  velocity  V  at  infinity. 

The  flow  about  a  circle  can  be  transformed  into  the  flow  around  an  aerofoil-like  shape  by 
substituting  the  variable 


-  32  - 


L 


into  the  expression  for  the  flow  about  a  circular  cylinder  having  a  radius  slightly  larger  than  a  (Rid.  2, 
p.  50).  The  necessary  placement  of  the  singularities  of  the  transformation  is  given  in  Reference  2.  This 
result  makes  the  general  transformation  in  Equation  (1)  specific  to  the  current  problem. 

The  expression  for  the  flow  about  a  circular  cylinder  of  radius  equal  to  (a  +  e)  in  the  z  plane 
with  the  flow  at  an  angle  a  to  the  x  axis  at  infinity  is  (Ref.  2,  p.  52): 


(a  +  e)-  il'  „  (z  +  e) 

w  =  V  (z  +  e)exp(-iao)  +  -  exp(ia,)  +  —  tn  -  exp(-ia, 


A  point  f  in  the  wing  plane  has  a  corresponding  point  z  in  the  cylinder  plane.  This  point 
can  be  obtained  by  using  Equation  (1-5): 

z  =  f/2  ±  xXf/2)"- a2  <!- 

where  the  sign  is  chosen  to  make  z  lie  outside  the  circle.  The  solution  is  best  performed  numerically 
because  f  and  z  are  complex. 

The  velocity  in  the  wing  plane  is  obtained  by  differentiating  Equation  (1-6)  with  respect 
to  z,  calculating  dz/df  from  Equation  (1-5)  and  substituting  the  results  in  Equation  (1-4): 


/  (a  +  e)2 

'  -  — r 


exp(io  ))  +  - 

0  /  27r(z  +  e) 


Note  that  this  result  is  in  terms  of  z.  Substitution  of  (1-7)  in  (1-8)  could  be  performed  but 
the  result  for  the  velocity  in  terms  of  f  is  unduly  cumbersome.  The  factor  outside  the  square  brackets 
comes  from  Equation  (1-5).  Equation  (1-8)  describes  the  flow  about  a  symmetrical  wing  whose  thick¬ 
ness  is  related  to  e. 

The  Kutta  Condition  (Ref.  3,  p.  68)  is  used  to  fix  the  rear  stagnation  point  to  the  wing 
trailing  edge  by  chosing  F  in  such  a  way  that  the  factor  inside  the  square  brackets  in  Equation  (1-8) 
is  equal  to  zero  at  z  =  a.  The  result  is 

T  =  47r(a  +  e)  Vsin  a  (1-9) 


APPENDIX  II 


LOCATION  OF  PROBE  SAMPLE  CENTERS 


Chord  =  78”  a  =  19.5  inches 

Distance  aft  of  leading  edge:  7.90  inches  £/a  =  -1.595 

Distance  below  chord  line:  short  pylon,  11.31  inches;  77/a 

long  pylon,  15.56  inches;  77, 'a 
Position  uncertainty  ±  3/16";  £/a  and  77/a  =  ±  .01 

Spanwise  location:  Outboard  s  =  .747 

Inboard  s  =  .618 
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SUMMARY/SOMMAIRE 

The  N  A I  Twin  Otter  atmospheric  research  aircraft  carries  instruments  mounted  under  the  wines  to  count  and  image  cloud  and 
precipitation  particles  for  cloud  physics  studies.  As  part  of  an  investigation  on  the  influence  of  the  aircraft  presence  on  these  measure¬ 
ments.  airflow  velocities  about  the  wing  have  been  calculated.  These  theoretical  results  have  been  compared  with  in-flight  measurements 
taken  by  pressure  probes  mounted  in  a  cannister  suspended  from  a  long  and  a  short  pylon  below  the  wing  at  two  spanwise  locations. 

The  Twin  Otter  wing  has  been  represented  by  a  Joukowski  aerofoil  of  the  same  maximum  thickness  but  w  it  It  a  cusped  trailing 
edge.  The  wing  camber  has  been  neglected  but  aspect  ratio  and  plantorm  effects  have  been  included.  The  elTcc*  of  the  cannister  was 
included  in  the  comparison  of  theory  with  experiment.  I' sc  of  the  theory  aft  of  the  quarter  chord  will  require  experimental  verification. 

Both  experimental  and  theoretical  results  for  tl.e  local  How  velocity  ratio  are  linear  functions  of  aircraft  lift  coefficient 
(CL>  in  the  range  of  Cl  between  0.35  to  1.06.  The  worst  discrepancy  between  theory  and  experiment  in  theCj  range  of  0.4  to  0.6  is 
-2.8'-  to  +1 .1'?.  The  results  for  the  inboard/long  pylon  and  outboard/short  pylon  have  considerably  less  deviation. 

The  theory  has  been  experimentally  verified  to  the  above  accuracy  and  can  be  used  to  describe  the  spatial  variation  ot  the 
How  velocity  ahead  of  the  wing.  Applications  of  the  theory  to  the  study  of  aircraft  How  effects  on  cloud  physics  measurements  can 
proceed  with  a  quantifiable  level  of  confidence. 
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